Quantum Radiation from Quantum Gravitational Collapse 
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We study quantum radiation emitted during the collapse of a quantized, gravitating, spherical do- 
main wall. The amount of radiation emitted during collapse now depends on the wavefunction of the 
collapsing wall and the background spacetime. If the wavefunction is initially in the form of a sharp 
wavepacket, the expectation value of the particle occupation number is determined as a function 
of time and frequency. The results are in good agreement with our earlier semiclassical analysis 
and show that the quantum radiation is non-thermal and evaporation accompanies gravitational 
collapse. 
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Quantum radiation from gravitationally collapsing 
matter has been the subject of a lot of interest (see e.g. 
[EBB 0,01 for early work and for more recent re- 

views) . Virtually all of the work in this context has been 
done in the semiclassical approximation in which the col- 
lapsing matter is described in classical terms and only the 
radiation field is quantized. Such results, though very im- 
portant, leave open the possibility that quantum collapse 
could qualitatively change the outcome. This would be 
relevant to both the end-point of Hawking radiation as 
well as to the dynamics of gravitational collapse close to 
the point of black hole formation. However, a full quan- 
tum treatment of gravitational collapse has been lacking 
so far. 

In Ref. 01 a functional Schrodinger formalism ap- 
proach to Hawking evaporation problem was developed, 
and the quantum evaporation from a classical collapsing 
background was studied. This formalism, being inher- 
ently quantum mechanical, is particularly convenient for 
treating the collapsing gravitational background in the 
context of quantum mechanics. In this paper we address 
the problem where the infalling matter and the radiation 
field are both treated in quantum theory. 

Our approach follows Ref. [§] and uses the functional 
Schrodinger equation 
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where "J is the wavefunctional for a collapsing spherical 
domain wall and the excitation modes of the radiation 
field. We will find the Hamiltonian more explicitly below, 
after which we will solve the equation for the wavefunc- 
tional. The wavefunctional is used to find the excitation 
spectrum of the radiation field which leads to the flux of 
quantum radiation from the collapsing matter. 

To study a concrete realization of black hole formation 
we consider a spherical Nambu-Goto domain wall that is 
collapsing. To include the possibility of (spherically sym- 
metric) radiation we consider a massless scalar field, 4>, 
that is coupled to the gravitational field but not directly 



to the domain wall. The action for the system is 



S = Id x^J —g 
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where the first term is the Einstein- Hilbcrt action for the 
gravitational field, the second is the scalar field action, 
the third is the domain wall action in terms of the wall 
world volume coordinates, £ a (a = 0, 1, 2), the wall ten- 
sion, c, and the induced world volume metric 



lab = gavdaX^dbX" 
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The coordinates X M (£ a ) describe the location of the wall, 
and Roman indices go over internal domain wall world- 
volume coordinates ( a , while Greek indices go over space- 
time coordinates. 

We will now only consider spherical domain walls and 
assume spherical symmetry for the scalar field ($ = 
<&(t, r)). Then the wall is described by only the radial de- 
gree of freedom, R(t). Furthermore, the metric is taken 
to be the solution of Einstein equations for a spherical 
domain wall. The metric is Schwarzschild outside the 
wall, as follows from spherical symmetry [lol ] 



ds 2 
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where, Rs = 2GM is the Schwarzschild radius in terms of 
the mass, M, of the wall, and d£l 2 is the usual angular line 
element. In the interior of the spherical domain wall, the 
line element is flat, as expected by Birkhoff's theorem, 



ds 2 = -dT 2 + dr 2 
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r < R(t) (5) 



The interior time coordinate, T, is related to the observer 
time coordinate, t, via the proper time, r, of the domain 
wall. 



dT _ 
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The ratio of these equations connects the interior and 
exterior time coordinates 



dT 
~dT 
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where R T = dR/dr and R = dR/dt. Integrating Eq. (JSj> 
still requires knowing R(t) (or R(t)) as a function of r 
(or t). 

Note that the wall radius, R(t), completely determines 
the metric. Hence, when we quantize R, both the wall 
and the background will be described quantum mechan- 
ically, but the wall and the background will not have 
independent quantum dynamics. 

Regarding our use of the Schwarzschild time coordi- 
nate, it is important to point out that the Schwarzschild 
line element has a coordinate singularity at r = Rs- 
However, the restriction in Eq. JJ to r > R(t) means 
that our coordinate system is well-defined everywhere as 
long as R(t) > Rs- In all that follows, we will only be 
considering the case R(t) > Rs- For example, we will 
calculate the quantum radiation from a collapsing wall 
with radius R(t) > Rs- 

The energy function for the wall is Q 



wall 
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where overdots denote derivatives with respect to t. This 
is a first integral of the equations of motion and is identi- 
fied with the energy of the gravitating domain wall fioj . 
The canonical momentum is given by 



n 
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where fx = cr(l — 2ttG<jRs)- Then, in the regime R • 
the wall Hamiltonian is 



H 
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and has the form of the energy of a relativistic particle, 



with a position dependent mass. In the limit 



B ~ 0, the mass term can be neglected 
ultra-relativistic - and hence 



H. 
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the wall is 



(13) 



where we have chosen the negative sign appropriate for 
describing a collapsing wall. 

The scalar field, $, is decomposed into a complete set 
of basis functions denoted by {fk{r)} 



k 
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The exact form of the functions fk(r) will not be impor- 
tant for us. We will be interested in the wavefunction for 
the mode coefficients {a^}. 

The Hamiltonian for the scalar field modes is found 
by inserting the scalar field mode decomposition and the 
background metric into the action 



= J d A x^^d^d^ 
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The Hamiltonian for the scalar field modes takes 
the form of coupled simple harmonic oscillators with 
R— dependent mass and couplings due to the non-trivial 
metric. In the regime R ~ i?Sj_for a normal mode de- 
noted by b, the Hamiltonian is 



Rs 



H h = 1 - — + —b 
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where it is the momentum conjugate to 6, m and K are 
constants whose precise values are not important for us. 

Hence the total quantum Hamiltonian for the wall and 
a normal mode of the scalar field is 



H = H, 
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where n = —id/dR, it = —id/db. The wavefunction is 
now a function of b, R and t: ^ = ^(b, R, t). The total 
Hamiltonian (|17[) is Hermitian with respect to the inner 
product with measure l/B. 

The differential form of the Schrodingcr equation is 
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dR 2m db 2 ' 2 
We look for stationary solutions 

f(b,R,t) = e- lEt tfj(b,R) (19) 

Then the time-independent Schrodinger equation is 
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where oj 2 = k/(mB) and e = E/B. To solve Eq. (|2"U| we 
use the ansatz 



i>(b,R) = e- i J R < R ' )dR 'ct>(b,R) = 
where 
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This leads to the equation for </>(&, 77) where r\ = Rs — R 
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As discussed in Ref. ll| , this has the implicit solution 

1/4 



4>(b,v) 
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where p r] is the derivative of the function p(j]) with re- 
spect to 77, and the defining equation for p is 



p m +uj 2 {rj)p = i 



where 
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The initial conditions for p are taken at some large neg- 
ative value of rj (i.e. large value of R) denoted by rji, so 
that 



p(m) 



1 
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The phase a is defined by 

1 f v drj 



a(r\) 



2 J p2(r)>) 
Then Eqs. HU) and (21]) give 
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with </> given in Eq. (|2~4"|) and, as above, r/ = R$ — R. 

Next we note that -E only enters the solution in the 
first exponential factor. So we can superpose stationary 
solutions to construct wavepackets. For example, 



*(b,R,t) = 



1 
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where a is the width of the wavcpackct, is a solution 
to the time dependent Schrodingcr problem in Eq. (jT5J). 
With the chosen initial conditions for p at rji, this solu- 
tion describes a wavepacket for the collapsing wall and 
quantum fields in their ground state at r\i. 

Note that we have normalized the wavepacket with 
unit measure in the u coordinate since the Schrodingcr 
evolution preserves this normalization. Also, u can be 
expressed in terms of r\ through Eq. (f2"2"f . 

The solution (|30p describes quantum radiation from a 
quantized collapsing shell of matter. The shell is repre- 
sented by a wave packet that has a Gaussian fall-off and 
is moving toward the horizon Rs which corresponds to 
u — > — oo. (The wave packet is not strictly Gaussian be- 
cause the function <fi also depends on R via rj.) Radiation 
from the shell, represented by <f>(b, 77) depends on the po- 
sition of the shell. We want to find the main features of 
this radiation. 

It is worth noting that 



*{b,R,t) = f(u + t)4>{b,r,) 



(31) 



is a solution for any function /, provided <p(b, 77) satisfies 
Eq. ((23)) . This solution shows that wavepackets in u do 
not spread in time at late times. Since dR = Bdu from 
Eq. (12"2")) and B gets smaller as the wall collapses, the 
wavepacket gets more sharply peaked in R with time. 

In the semiclassical analysis where R is treated clas- 
sically, the wavefunction for the mode amplitudes is de- 
composed into suitably chosen basis wavefunctions (dis- 
cussed below). If a complete basis is denoted by {0 n }, 
then the expectation of the occupation number is 



(32) 



Similar considerations hold for the present case where 
we are treating R in quantum theory. The only difference 
is that the value of N depends on what value R takes. 
Hence N is given by an R— and t— dependent probability 
distribution, and to find the expectation value of N at 
any time, we also need to integrate over R. Therefore, 



du (N(R, t)) = / duJ2n\(<j> n \y)\ 2 (33) 



(N(t)) 



Here we have been careful to integrate over du with unit 
measure which is equivalent to integrating over R but 
with measure 1/B. 

To proceed further, we need to specify our basis wave- 
functions, 0„. These are chosen to be simple harmonic 
oscillator basis states at a frequency u) 



Mb) 



fmu> 

\ 7T 



-mab 2 /2 



V2"nl 

where H n are Hermite polynomials. Then Q 
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(34) 



6„|*) =f(u+t)- 



If lY /2 ("-!)'! 
P ~ P J ~~J^~ 



(35) 

where f(u + t) can be chosen to be a Gaussian function 
as in Eq. (201) and 



P= 1- 



Pv 



(36) 



The sum over n in Eq. 
and the occupation number is 



VP p- 

can be done explicitly § 



(N(t)) 



du\f(u + t)\ 
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Note that the integrand contains p which is a function 
of 77 = Rs — R (Eq. I|25p). which in turn is a function 
of u (Eq. ([22]) ). The frequency uj in Eq. |26|) depends 
on R and hence is a quantum variable. The frequency Q 
at which N is evaluated is chosen to be the expectation 
value of to: Q — (uj). 
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FIG. 1: ln(JV) versus t/Rs for various fixed values of ujRs 
and for the wavepacket width a = 0.1i?s- Here, uj is the 
expectation value of the frequency calculated for the state 
(|30l) . The linear growth shows that the occupation number 
increases exponentially at late times. 
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t/R s = ll, 12,13 




FIG. 2: N versus uiRs for various fixed values of t/Rs for 
a = 0.1-Rs- The oscillations depend on the choice of width 
of the wavepacket a. Here, u> is the expectation value of 
the frequency calculated for the state (|30[) . The occupation 
number at any frequency grows as t/Rs increases. 



To evaluate (TV) we now need to solve Eq. (|25p and then 
insert the solution in the integrand in Eq. (|37[) together 
with the Gaussian form of f(u+t). Finally, we need to do 
the u integral in Eq. l[57|). If f(u + t) is sharply peaked, 
only a small window of u, and an even smaller window 
of R (Eq. contributes to the integral. 

We have calculated the occupation number numerically 
for the case of a wavepacket with a = O.IRs and we show 
the plot of (N(t)) versus t in Fig. [TJ These calculations, 
when done for a range of (D, also yield the spectrum of 
occupation numbers. In Fig. [5] we show the spectrum for 
a few different times. This result has the same general 
features as in the semiclassical analysis |{|. The main 
effect of the quantized background is that the spectrum 
becomes more non-thermal though this depends on the 
width of the wavepacket. In the case of a very sharp 
wavepacket (<x — ► 0), the spectrum becomes identical to 
the semiclassical spectrum. 

In summary, we have studied quantum radiation from 
the changing spacetime metric due to a collapsing spheri- 
cal domain wall. The main improvement over earlier cal- 



culations is that the collapsing matter i.e. domain wall 
is treated within quantum theory. The end results for 
the growth of particle occupation number and spectrum 
are very similar to that found in semiclassical analysis 
and, in fact, become identical to the semiclassical result 
as the width of the wavepacket decreases. 

Our analysis implies that quantum radiation is emitted 
continuously during quantum gravitational collapse, and 
so evaporation and collapse are concomitant. As seen 
in Fig. [TJ the occupation number at any frequency grows 
with time, and as in Fig. [21 the spectrum has non-thermal 
features. In an earlier paper we have approximately 
fitted a thermal distribution to the spectrum in an in- 
termediate range of frequencies, and found a radiation 
roughly consistent with the Hawking temperature. As 
the collapse proceeds, the window of frequencies in which 
the radiation is thermal grows, but only in the t — > oo 
limit, does the spectrum become thermal in an infinite 
range of frequencies. The thermal spectrum, however, 
may never be achieved in a physical setting precisely be- 
cause it is only realized in the t — > oo limit. It is also very 
important that small non-thermalities found in the semi- 
classical treatment Q get amplified in the full quantum 
treatment. 

We should note that our analysis has ignored radiation 
backreaction on the collapse process since we take Rs to 
be constant. A full study of backreaction must involve 
time dependence of Rs- 
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